\(\int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx\) [150]

Optimal result

Integrand size = 25, antiderivative size = 66 \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\frac {x^2 (1-a x)}{a^2 \sqrt {1-a^2 x^2}}+\frac {(4-3 a x) \sqrt {1-a^2 x^2}}{2 a^4}+\frac {3 \arcsin (a x)}{2 a^4} \]

[Out]

3/2*arcsin(a*x)/a^4+x^2*(-a*x+1)/a^2/(-a^2*x^2+1)^(1/2)+1/2*(-3*a*x+4)*(-a^2*x^2+1)^(1/2)/a^4

Rubi [A] (verified)

Time = 0.03 (sec) , antiderivative size = 66, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.160, Rules used = {864, 833, 794, 222} \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\frac {3 \arcsin (a x)}{2 a^4}+\frac {x^2 (1-a x)}{a^2 \sqrt {1-a^2 x^2}}+\frac {(4-3 a x) \sqrt {1-a^2 x^2}}{2 a^4} \]

[In]

Int[x^3/((1 + a*x)*Sqrt[1 - a^2*x^2]),x]

[Out]

(x^2*(1 - a*x))/(a^2*Sqrt[1 - a^2*x^2]) + ((4 - 3*a*x)*Sqrt[1 - a^2*x^2])/(2*a^4) + (3*ArcSin[a*x])/(2*a^4)

Rule 222

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Simp[ArcSin[Rt[-b, 2]*(x/Sqrt[a])]/Rt[-b, 2], x] /; FreeQ[{a, b}
, x] && GtQ[a, 0] && NegQ[b]

Rule 794

Int[((d_.) + (e_.)*(x_))*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[((e*f + d*g)*(2*p
+ 3) + 2*e*g*(p + 1)*x)*((a + c*x^2)^(p + 1)/(2*c*(p + 1)*(2*p + 3))), x] - Dist[(a*e*g - c*d*f*(2*p + 3))/(c*
(2*p + 3)), Int[(a + c*x^2)^p, x], x] /; FreeQ[{a, c, d, e, f, g, p}, x] &&  !LeQ[p, -1]

Rule 833

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Simp[(d + e*x)^(m
- 1)*(a + c*x^2)^(p + 1)*((a*(e*f + d*g) - (c*d*f - a*e*g)*x)/(2*a*c*(p + 1))), x] - Dist[1/(2*a*c*(p + 1)),
Int[(d + e*x)^(m - 2)*(a + c*x^2)^(p + 1)*Simp[a*e*(e*f*(m - 1) + d*g*m) - c*d^2*f*(2*p + 3) + e*(a*e*g*m - c*
d*f*(m + 2*p + 2))*x, x], x], x] /; FreeQ[{a, c, d, e, f, g}, x] && NeQ[c*d^2 + a*e^2, 0] && LtQ[p, -1] && GtQ
[m, 1] && (EqQ[d, 0] || (EqQ[m, 2] && EqQ[p, -3] && RationalQ[a, c, d, e, f, g]) ||  !ILtQ[m + 2*p + 3, 0])

Rule 864

Int[((x_)^(n_.)*((a_) + (c_.)*(x_)^2)^(p_))/((d_) + (e_.)*(x_)), x_Symbol] :> Int[x^n*(a/d + c*(x/e))*(a + c*x
^2)^(p - 1), x] /; FreeQ[{a, c, d, e, n, p}, x] && EqQ[c*d^2 + a*e^2, 0] &&  !IntegerQ[p] && ( !IntegerQ[n] ||
!IntegerQ[2*p] || IGtQ[n, 2] || (GtQ[p, 0] && NeQ[n, 2]))

Rubi steps \begin{align*} \text {integral}& = \int \frac {x^3 (1-a x)}{\left (1-a^2 x^2\right )^{3/2}} \, dx \\ & = \frac {x^2 (1-a x)}{a^2 \sqrt {1-a^2 x^2}}-\frac {\int \frac {x (2-3 a x)}{\sqrt {1-a^2 x^2}} \, dx}{a^2} \\ & = \frac {x^2 (1-a x)}{a^2 \sqrt {1-a^2 x^2}}+\frac {(4-3 a x) \sqrt {1-a^2 x^2}}{2 a^4}+\frac {3 \int \frac {1}{\sqrt {1-a^2 x^2}} \, dx}{2 a^3} \\ & = \frac {x^2 (1-a x)}{a^2 \sqrt {1-a^2 x^2}}+\frac {(4-3 a x) \sqrt {1-a^2 x^2}}{2 a^4}+\frac {3 \sin ^{-1}(a x)}{2 a^4} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.19 (sec) , antiderivative size = 69, normalized size of antiderivative = 1.05 \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\frac {\sqrt {1-a^2 x^2} \left (4+a x-a^2 x^2\right )}{2 a^4 (1+a x)}+\frac {3 \arctan \left (\frac {a x}{-1+\sqrt {1-a^2 x^2}}\right )}{a^4} \]

[In]

Integrate[x^3/((1 + a*x)*Sqrt[1 - a^2*x^2]),x]

[Out]

(Sqrt[1 - a^2*x^2]*(4 + a*x - a^2*x^2))/(2*a^4*(1 + a*x)) + (3*ArcTan[(a*x)/(-1 + Sqrt[1 - a^2*x^2])])/a^4

Maple [A] (verified)

Time = 0.37 (sec) , antiderivative size = 97, normalized size of antiderivative = 1.47

method	result	size
risch	\(\frac {\left (a x -2\right ) \left (a^{2} x^{2}-1\right )}{2 a^{4} \sqrt {-a^{2} x^{2}+1}}+\frac {3 \arctan \left (\frac {\sqrt {a^{2}}\, x}{\sqrt {-a^{2} x^{2}+1}}\right )}{2 a^{3} \sqrt {a^{2}}}+\frac {\sqrt {-\left (x +\frac {1}{a}\right )^{2} a^{2}+2 \left (x +\frac {1}{a}\right ) a}}{a^{5} \left (x +\frac {1}{a}\right )}\)	\(97\)
default	\(\frac {\arctan \left (\frac {\sqrt {a^{2}}\, x}{\sqrt {-a^{2} x^{2}+1}}\right )}{a^{3} \sqrt {a^{2}}}+\frac {-\frac {x \sqrt {-a^{2} x^{2}+1}}{2 a^{2}}+\frac {\arctan \left (\frac {\sqrt {a^{2}}\, x}{\sqrt {-a^{2} x^{2}+1}}\right )}{2 a^{2} \sqrt {a^{2}}}}{a}+\frac {\sqrt {-a^{2} x^{2}+1}}{a^{4}}+\frac {\sqrt {-\left (x +\frac {1}{a}\right )^{2} a^{2}+2 \left (x +\frac {1}{a}\right ) a}}{a^{5} \left (x +\frac {1}{a}\right )}\)	\(134\)

[In]

int(x^3/(a*x+1)/(-a^2*x^2+1)^(1/2),x,method=_RETURNVERBOSE)

[Out]

1/2*(a*x-2)*(a^2*x^2-1)/a^4/(-a^2*x^2+1)^(1/2)+3/2/a^3/(a^2)^(1/2)*arctan((a^2)^(1/2)*x/(-a^2*x^2+1)^(1/2))+1/
a^5/(x+1/a)*(-(x+1/a)^2*a^2+2*(x+1/a)*a)^(1/2)

Fricas [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 75, normalized size of antiderivative = 1.14 \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\frac {4 \, a x - 6 \, {\left (a x + 1\right )} \arctan \left (\frac {\sqrt {-a^{2} x^{2} + 1} - 1}{a x}\right ) - {\left (a^{2} x^{2} - a x - 4\right )} \sqrt {-a^{2} x^{2} + 1} + 4}{2 \, {\left (a^{5} x + a^{4}\right )}} \]

[In]

integrate(x^3/(a*x+1)/(-a^2*x^2+1)^(1/2),x, algorithm="fricas")

[Out]

1/2*(4*a*x - 6*(a*x + 1)*arctan((sqrt(-a^2*x^2 + 1) - 1)/(a*x)) - (a^2*x^2 - a*x - 4)*sqrt(-a^2*x^2 + 1) + 4)/
(a^5*x + a^4)

Sympy [F]

\[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\int \frac {x^{3}}{\sqrt {- \left (a x - 1\right ) \left (a x + 1\right )} \left (a x + 1\right )}\, dx \]

[In]

integrate(x**3/(a*x+1)/(-a**2*x**2+1)**(1/2),x)

[Out]

Integral(x**3/(sqrt(-(a*x - 1)*(a*x + 1))*(a*x + 1)), x)

Maxima [A] (verification not implemented)

none

Time = 0.29 (sec) , antiderivative size = 68, normalized size of antiderivative = 1.03 \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\frac {\sqrt {-a^{2} x^{2} + 1}}{a^{5} x + a^{4}} - \frac {\sqrt {-a^{2} x^{2} + 1} x}{2 \, a^{3}} + \frac {3 \, \arcsin \left (a x\right )}{2 \, a^{4}} + \frac {\sqrt {-a^{2} x^{2} + 1}}{a^{4}} \]

[In]

integrate(x^3/(a*x+1)/(-a^2*x^2+1)^(1/2),x, algorithm="maxima")

[Out]

sqrt(-a^2*x^2 + 1)/(a^5*x + a^4) - 1/2*sqrt(-a^2*x^2 + 1)*x/a^3 + 3/2*arcsin(a*x)/a^4 + sqrt(-a^2*x^2 + 1)/a^4

Giac [F(-2)]

Exception generated. \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate(x^3/(a*x+1)/(-a^2*x^2+1)^(1/2),x, algorithm="giac")

[Out]

Exception raised: TypeError >> an error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:sym2poly/r2sym(const gen & e,const index_m & i,const vecteur & l) Error: Bad Argument Value

Mupad [B] (verification not implemented)

Time = 0.07 (sec) , antiderivative size = 116, normalized size of antiderivative = 1.76 \[ \int \frac {x^3}{(1+a x) \sqrt {1-a^2 x^2}} \, dx=\frac {3\,\mathrm {asinh}\left (x\,\sqrt {-a^2}\right )}{2\,a^3\,\sqrt {-a^2}}-\frac {\left (\frac {1}{a^2\,\sqrt {-a^2}}+\frac {x\,\sqrt {-a^2}}{2\,a^3}\right )\,\sqrt {1-a^2\,x^2}}{\sqrt {-a^2}}-\frac {\sqrt {1-a^2\,x^2}}{a^3\,\left (x\,\sqrt {-a^2}+\frac {\sqrt {-a^2}}{a}\right )\,\sqrt {-a^2}} \]

[In]

int(x^3/((1 - a^2*x^2)^(1/2)*(a*x + 1)),x)

[Out]

(3*asinh(x*(-a^2)^(1/2)))/(2*a^3*(-a^2)^(1/2)) - ((1/(a^2*(-a^2)^(1/2)) + (x*(-a^2)^(1/2))/(2*a^3))*(1 - a^2*x
^2)^(1/2))/(-a^2)^(1/2) - (1 - a^2*x^2)^(1/2)/(a^3*(x*(-a^2)^(1/2) + (-a^2)^(1/2)/a)*(-a^2)^(1/2))